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The formation of patterned nanostructure in AgCo/Pt/MgO(001) thin film is simulated by a 
technique of combining molecular dynamics and phase-field theory. The dislocation (strain) network 
existing in Pt/MgO is used as a template whose pattern is transferred to AgCo phase in spinodal 
decomposition, resulting in regular arrays of Co islands that are attracted by the dislocations. 
The influence of various factors, such as component concentration and film thickness, is studied, 
ft is found that the spinodal decomposition of AgCo in this system is mainly characterized by a 
competition between a surface-directed layer structure and the strain-induced patterned structure, 
where the patterned Ag-Co structure only dominates in a small range near the interface (less than 
10 atomic layers). However, if the interlayer diffusion can be minimized by controlling film growth 
conditions, it is shown that the patterned structure can be formed throughout the entire film. 

PACS numbers: 68.55.-a, 64.75.+g, 81.16.Rf, 81.80.Hd 



I. INTRODUCTION 

With the rapid development of the high-density infor- 
mation storage device, the control of grain size and grain 
size distribution in magnetic thin films will be required 
for the next generation of information storage technol- 
ogy, where the drive for decreased media noise and in- 
creased storage density is pushing the grain size below 
the 10 nm regime. [Q p| Present technology involves the 
deposition of Co rich magnetic phases onto seed layers 
whose grain size determines the grain size of the mag- 
netic layer through grain to grain epitaxy. For example, 
in IBM's 10 Gbit/in. 2 demonstration, the magnetic film 
with an average grain size of 12 nm is obtained where 
a Co phase is deposited onto NiAl seed layers. Q Other 
methods to control grain size and grain size distribution 
were also explored in the labs. It was shown recently that 
monodispersed magnetic nanoparticles can be created by 
chemical means and deposited onto a substrate to form 
a highly regular array of particles. H 

There are several possibilities for introducing control- 
lable length-scales into the film formation process, one 
of which is the self-organized growth on a substrate 
with preferred nano-pattern. For example, flat nano- 
sized Co dots form in Au(lll) substrate with surface 
reconstruction, Q and ordered arrays of vertical magnetic 
Co pillars can be further obtained by alternative deposi- 
tion of Co and Au layers. (j| When Co film is deposited on 
the Pt(lll), fee and hep phases are arranged in regular 
patterns to reduce the strain energy. Q Recently, Kern 
et al showed that the strain-relief dislocation network in 
substrate can be transferred through heterogeneous nu- 



cleation to regular superlattice of almost mono-dispersed 
islands with great feasibility. 

In a previous paper, [|| the interaction between the 
spinodal decomposition of thin film and the disloca- 
tion network of substrate was briefly discussed with 
AgCo/Pt/MgO thin film as an example. It was re- 
vealed that Co phases are attracted by dislocations in 
the Pt/MgO interface, producing an interesting nanos- 
tructure. Such mechanism may provide a possibility to 
control the length scale for magnetic recording media 
layer. For example, after the patterned nanostructure 
of AgCo phases is established in spinodal decomposition 
as a seed layer, a Co-rich conventional medial alloy can 
grow on it. In this article, detailed simulations are con- 
ducted in AgCo/Pt/MgO(001) thin film to investigate 
the formation of nanostructure under various conditions. 
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FIG. 1: Schematic graphics of the AgCo/Pt/MgO(001) thin 
film. The symbols cross (+), triangle (A) and stars (*) rep- 
resent atoms Mg, O, and Pt, respectively, while filled circle 
(•) for AgCo atoms. Pt atoms trend to sit atop the sur- 
face O atoms due to their favorable interaction energy. The 
arrows indicate the locations where dislocation will be gener- 
ated. Axes are plotted in units of A. 
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FIG. 2: The strain effect of the as-grown Pt/MgO(001) thin 
film that is represented by non-uniform atomic displacements 
on the surface. A system with 3 ML of Pt is simulated and 
the displacements of surface atoms are plotted. 



cle, we adopt an alternative method ||, [L3| of combining 
atomistic calculations and microscopic mean-field theory, 
which is applicable for noncoherent cases, to simulate the 
effect of dislocation, strain response and spinodal decom- 
position in AgCo/Pt/MgO(001) thin film. In this study, 
the morphology of the system is described by the position 
and occupation probability (atomic concentration) of all 
atoms, {i"i,ni}. {r^} are used to simulate local atomic 
motion such as strain response or dislocation, while {rii} 
are used to simulate atomic diffusion of phase transforma- 
tion. The atomic concentrations of MgO and Pt phases 
are fixed in simulation, while the AgCo phase undergoes 
a process of spinodal decomposition. After expressing 
the empirical potential and free energy as functions of 
{r^, rii} under mean-field approximation, i.e., V({ri, rij}) 
and _F({r,-, n^}), the evolution of atomic positions is given 
in molecular dynamics (MD) asflhft 



II. THEORETICAL MODEL AND SIMULATION 
METHOD 

A. System model 

The system considered here is the AgCo/Pt/MgO(00I) 
thin film (see Fig. 1). MgO is preferred substrate material 
in many of the new information storage and processing 
devices made from thin metal films. || MgO has a f.c.c. 
lattice structure with lattice parameter aM g o = 

4.2IA, 

which is slightly larger than that of Pt (apt = 3.924A). 
When Pt film is deposited on MgO substrate, due to the 
interactions between film and substrate, Pt phase will 
be subject to elastic mismatch strain. The mismatch 
strain energy is partly released by the formation of inter- 
facial dislocation network. jl(| The dislocation network of 
Pt/MgO(00I) with cube-on-cube orientation relationship 
appears as a square net with dislocation spacing 4.05 nm. 
Such dislocation network results in periodic variations of 
strain in Pt thin film (see Fig. 2). ft has been pointed out 
by theoretical analysis|TlJ and simulation^ that such pe- 
riodic strain provides a nanopattern for the spinodal of 
a second film on it. fn this article, we consider a AgCo 
thin film on the strained Pt/MgO(001) (Fig. I) and sim- 
ulate the spinodal decomposition of AgCo phase. As 
most simulations on spinodal decomposition, the atoms 
of the initial AgCo phase are randomly arranged in regu- 
lar (f.c.c.) crystal lattice. The model assumes that AgCo 
is deposited on Pt/MgO in a disordered state and we cal- 
culate the formation of nanostructure by a subsequent 
annealing process. 



B. Simulation method 

There are a number of approaches proposed to model 
the elastic effect on precipitate morphology in coher- 
ent systems, jl2| where the linear elasticity theory is 
employed to calculate the elastic energy. In this arti- 



™^ = -v^({iW), (f) 



while the diffusion pro cess is described in microscopic 
mean- field theory asjf5|, |f6[| 

3 J 

where Lij is the kinetic coefficient proportional to the in- 
verse average time of elementary diffusional jumps from 
site i to j. The timescale of diffusion is mainly deter- 
mined by the magnitude of L^. In our calculation, L t j 
is assigned as a constant L when — r j- 1 < i?diff and 
otherwise. Since diffusion is much slower than elastic 
response, the system is set in elastic equilibrium during 
the diffusion process. In numerical scheme, the simula- 
tion loop is as following: 

(1) run m steps of MD by Eq. (I) with time step St to 
reach elastic equilibrium; 

(2) run one step of diffusion by Eq. (2) with time step dt; 

(3) goto (I). 

It is a multiscale scheme overcoming the weakness of 
small time scale in MD since mSt <C dt. Although 
the present work employed a microscopic diffusion equa- 
tion for a binary system and assumed an atomic ex- 
change mechanism for diffusion, it is rather straightfor- 
ward to formulate the microscopic diffusion using a va- 
cancy mechanism by introducing vacancy as the third 
component. With vacancy as the third component, the 
model can automatically model the surface morphology 
changes as a result of surface energy anisotropy and the 
balance of surface energies and the interphase boundary 
energy between the two phases in the film. [r^| However, 
as the first attempt, we assumed the simple exchange 
mechanism. The vacancy mechanism will be considered 
in a future publication. 
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C. Empirical potential 

The empirical potential between Ag, Co and Pt 
atoms adopted in current simulation is the Tight-binding 
second- moment-approximation (TB-SMA) scheme, [18 
where the potential is comprised of a binding part and a 
repulsive part: 



For AgCo phase, to include the influence of atomic con- 
centration, the following formula is used under mean-field 
approximation for the potential (A and B denotes two 
atomic species, Ag and Co, respectively, and rii, nj are 
the concentrations of A atoms): 



(3) 



rijn 



j 



Aaa exp[-pAA(»»jAo _ !)] 
+ [m(l - nj) + (1 - nijrij] Aab exp[-pAB(rij /r^ B - 1)] 



+ (1 - n*) (1 - nj)A B B exp[-p BB {nj/r BB 



(4) 



1)] 



E B = 



E 

3 



niUj^expl-qAA^ij /rfr A - 1)] 
+ [rvi(l -rij) + (1 - TH)nj] £ AB ex P[ 
+ (1 - rii){l - nj)£ BB exp[-qBB(nj/r BB - 1)] 



1/2 



-g AB (r y /r AB - 1)] 



(5) 



where r»j represents the distance between atoms i and 
j and Tq' 3 is the first-neighbors distance in the a/3 lat- 
tice. A, p, £, q are some free parameters of the SMA 
scheme that are fitted to desired properties of the sys- 
tem. The parameters for the pure species are available 
in the literature. |Q For the cross parameters of Ag-Co 
interaction, we fit A, £ to the heat of formation (+26 
kJ/g atom)|l9) and lattice parameter (3.81 A that is as- 
sumed to be the average value between pure Ag and Co) 
of the Ago.sCoo.s disordered phase, while p and q are sim- 
ply assigned to the average values of pure species. The 
result obtained is shown in Table I. Based on the poten- 
tial scheme Eqns. (3-5) and the fitted parameters, the 
formation heat and the lattice constant at any atomic 
concentration can be calculated, which are depicted as 
Fig. 3. The curve of formation heat can be described by 
the parabolic law very well. The lattice constant approx- 
imately obeys the Vegard's law. The elastic constants 
are also calculated (see Fig. 4). It appears that the elas- 
tic constants monotonously change with concentration. 
However, it should be noted that the curves do not sat- 
isfy the linear law that is assumed by many spinodal de- 
composition simulations. 

The cross parameters of Ag-Pt and Co-Pt are also 
listed in Table I, which are fitted in a similar process. 
The heat of formation is assumed to be zero. 

The importance of the Pt/MgO substrate in the cur- 
rent simulation lies in its dislocation network and the 
corresponding strain effect. We keep the MgO atoms 
fixed in simulation and describe the interactions between 
Pt and MgO as the Lenard- Jones potential: 



The parameters are adopted as: ffMg-Pt = 0.336eV, 

fMg-Pt = 

2.3A, 

£o-Pt = 1.34eV, fo-pt = 

1.783 A, with 

a distance cutoff 3.4A. With these parameters, isolated 
Pt atoms can be adsorbed by MgO substrate by sitting 
1.99A atop the surface O atoms of MgO(OOl) with an 
adsorption energy of 2.35 eV as indicated by density- 
functional calculations. [^0| MD simulation shows that 
such potential can produce dislocation network as re- 
quired and the corresponding periodic strain (see Fig. 2). 
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FIG. 3: Heat of formation and lattice constant of disordered 
AgCo phase as functions of Ag concentration. The curves 
are calculated by using the mean-field TB-SMA potential 
[Eqns. (3-5)] with parameters in Table I. 
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TABLE I: TB-SMA parameters for Ag-Co-Pt 





A (eV) 


£ (eV) 


P 


q 


ro (A) 


Ag-Co 


0.0752 


1.1944 


11.261 


2.7127 


2.696 


Ag-Pt 


0.1988 


1.9399 


10.770 


3.5711 


2.832 


Co-Pt 


0.1449 


1.9959 


11.108 


3.145 


2.667 



III. RESULTS AND DISCUSSIONS 
A. Basic feature 

A 3D system of AgCo/Pt/MgO(001) thin film with 
1 monolayer (ML) of Pt and 10 ML of AgCo phase is 
simulated. Every monolayer (in xy-plane) contains 30 x 
30 AgCo atoms and periodic boundary conditions are 
applied in x- and y-axes. Free boundary condition is used 
in z-axis. The average Ag concentration is CA g = 0.5. 
The diffusion cutoff in one diffusion step is i?diff = 3. 53 A, 
i.e., Lij = L for the nearest neighbors and = in 
other cases. A reduced time is defined as t* = t/L to 
measure the simulation process. The time step of MD is 
kept as St — 0.15 x 10 _15 s, and the temperature is 600 
K. The time step of diffusion simulation is dynamically 
adjusted to ensure the convergence of the result. 

The simulation result corresponding to t* = 70 is de- 
picted in Fig. 5. 240000 MD steps are conducted in 
simulation. We can see that the surface of thin film 
is covered by Ag phase. The reason is that the sur- 
face energy of Ag (1.3 J/m 2 ) is half that of Co (2.71 
J/m 2 ).@ Similar case occurs in CuCo system where Cu 
phase will diffuse onto the surface to form a metastable 
alloy in epitaxial growth. ||^| Due to the influence of sur- 
face, AgCo decomposed into layer structure near the 
surface. In fact, it is a kind of surface-directed spin- 




0.2 0.4 0.6 0.8 

Ag concentration 

FIG. 4: Elastic constants of disordered AgCo phase as func- 
tions of Ag concentration under the mean-field TB-SMA po- 
tential. 




FIG. 5: Atomic configuration of the system with 1 ML of Pt 
and 10 ML of AgCo (c Ag = 0.5) at t* = 70. 240000 molecular 
dynamics steps are used. The symbols +, A and * represent 
Mg, O, and Pt atoms, respectively. The occupation probabil- 
ity rii of AgCo atoms is visualized by the filling darkness of 
the circle at n, with white representing low value (Co) and 
black representing high value (Ag). (a) The x-z monolayer 
located at y — 0. (b) The third x-y monolayer of AgCo phase 
(counting from the bottom to the surface). 



odal decomposition. |23|] Interesting phenomenon is ob- 
served near the AgCo/Pt interface: Co atoms are at- 
tracted to the strained regions near the dislocations be- 
cause their atomic radius are smaller than Ag. As a re- 
sult, AgCo decomposed phase produces interesting 2D 
nanostructure pattern [Fig. 5(b)]. So the basic feature of 
spinodal decomposition in AgCo/Pt/MgO(001) system 
can be described as a competition between the surface- 
directed layer structure and the substrate-dislocation- 
induced nanopattern. These features are consistent with 
previous results in 2D approximation. M 



B. Effect of AgCo concentration 

In this section, we investigate the effect of AgCo con- 
centration on the spinodal decomposition. Two Ag con- 
centrations, CAg = 0.25 and 0.75, are considered, and the 
results are compared with that of CAg = 0.5 in Section 
III.A. 

A result of CA g = 0.25 is given in Fig. 6. The 
layer structure near the surface is similar to the case of 
CAg = 0.5. The difference occurs in the nanostructure 
near the interface: due to the reduction of Ag concentra- 
tion, Ag-rich regions near the interface shrink in volume, 
and the nano-pattern changes from Co islands in Ag ma- 
trix [Fig. 5(b)] into Ag islands in Co matrix [Fig. 6(b)]. 
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(a) 



FIG. 6: Snap-shot of phase morphologies of the system with 
1 ML of Pt and 10 ML of AgCo (c Ag = 0.25). (a) The x-z 
monolayer located at y — 0. (b) The second x-y monolayer of 
AgCo phase. 





(b) 
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FIG. 7: Atomic configuration of the system with 1 ML of 
Pt and 10 ML of AgCo (c Ag = 0.75). (a) The x-z monolayer 
located at y — 0. (b) The third x-y monolayer of AgCo phase. 
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FIG. 8: Atomic configuration of the system with 1 ML of 
Pt and 30 ML of AgCo (c Ag = 0.5). (a) The x-z monolayer 
located at y — 0. (b) The third x-y monolayer of AgCo phase. 



All Co atoms condense as cone-like islands on the dis- 
locations of the substrate. They are covered by Ag-rich 
phase that spread to the surface. When ca s = 0.5, the 
Co phases near the interface also assume the shape of 
cone-like islands [Fig. 5(a)], but the cones stand upside 
down, merging with the layer structure near the surface. 

It is noted that the Ag phase is important in forming 
the Co islands in the above process. If there is no Ag 
phase, Co thin film will grow in layer structure at small 
thickness. Ag phase acts as equivalent surfactant in as- 
sisting Co to decompose into islands. The present study 
suggests it may be possible to produce nanostructure by 
the following process. First, Co islands can be formed 
via spinodal decomposition of AgCo on periodic strained 
substrate. Subsequently, the Ag phase can be removed 
by chemical etching, leaving a regular array of Co islands. 
Experiments will be needed to test this conjecture. 



When the system is evolved in longer time, it is observed 
in simulation that the flat Ag islands get smaller and 
smaller. It means the surface Ag phase tends to absorb 
the Ag islands near the interface. 

Now, we come to the case of CA g = 0.75 (Fig. 7). Com- 
pared with the result of CA g = 0.5, there is no layer 
structure of Co phase near the surface in present case. 



C. Effect of AgCo phase thickness 

To explore the effect of AgCo phase thickness, a system 
with 30 ML of AgCo with CA g = 0.5 is simulated with 
the result in Fig. 8. The basic characteristics is similar 
to the above results: a competition between the layer 
structure near the surface and the patterned structure 
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FIG. 9: Pattern correlation function (solid line and circles) 
and layer Ag concentration (dashed line and triangles) of 
AgCo phase in the system with 1 ML of Pt and 30 ML of AgCo 
(cA g = 0.5). The pattern correlation function is normalized 
to make sure the self-correlation is equal to 1 (R mm = 1 for 
m — 3 in this case). 



near the interface. The patterned structure induced by 
the dislocations of the substrate decays within a distance 
of about ten monolayers while the layer structure presides 
in most regions. Surface energy is on the same order as 
the atomic cohesive energy while strain energy is higher 
order variation of the cohesive energy. So the surface 
effect dominates at the current film thickness. 

In order to quantitatively analyze the competition be- 
tween two kinds of structures, we define the pattern cor- 
relation function of two layers as: 



V • V 

v m v n 



Vr, 



IV., 



— R° 



(7) 



where V TO , V„ are vectors composed of atomic concen- 
trations within layers m and n. i?^„ is the calculation 
result when there is no correlation between two patterns: 



R° 



|V m I ■ |V„ 



(8) 



where V mfJl are components of V m and Mayer is the num- 
ber of components. 

The third AgCo layer has a representative pattern in- 
duced by the dislocation of substrate. We calculate the 
correction function of all AgCo layers with respect to the 
third one and plot the result in Fig. 9, together with 
the layer concentration. It is clearly demonstrated that 
surface-directed spinodal decomposition (dashed line) 
penetrates throughout the thin film while the pattern 
structure (solid line) only exists in the range of ten mono- 
layers near the interface. It suggests that the film thick- 
ness should be limited to small values if one is interested 
in producing nanostructure by spinodal decomposition 
on strained substrate. For example, when the thickness 
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FIG. 10: Atomic configuration of the system with 1 ML of 
Pt and 5 ML of AgCo (cA g = 0.5). (a) The x-z monolayer 
located at y — 0. (b) The third x-y monolayer of AgCo phase. 





x(A) 



FIG. 11: Atomic configuration of the system with 3 ML of Pt 
and 30 ML of AgCo (cAg = 0.5). The third x-y monolayer of 
AgCo phase is depicted here. 



of AgCo phase decreases to 5 ML, a regular array of Co 
islands can be obtained in simulation (see Fig. 10). 



D. Effect of Pt thickness 

The strain induced by the dislocation declines in Pt 
phase. If the thickness of Pt phase is large, the strain 
acting on AgCo phase will be too weak to induce nano- 
pattern. 

Fig. 11 displays the simulation results on the system 
with 3 ML of Pt and 30 ML of AgCo. It can be seen 
that there are only two Co islands within the simulation 
domain, in contrast to four islands in the previous cases. 
The other two Co islands, which locate at y — 30 A in 
the system with 1 ML of Pt, now disappear. It means 
the strain spreading through 3 ML of Pt is losing its 
capability of attracting Co atoms to form islands. 3 ML 
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(a) 
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x(A) 

FIG. 12: Atomic configuration of the thin film grown layer 
by layer where atomic diffusion is restricted within the same 
layers. The system has 1 ML of Pt and 10 ML of AgCo 
(cA g = 0.75). (a) The x-z monolayer located at y = 0. (b) 
The third x-y monolayer of AgCo phase. 

is the critical thickness of Pt phase to produce patterned 
structure in the current case. 

It has been illuminated in the previous sections that 
the influence of surface is competing with that of disloca- 
tion strain. With decreasing AgCo thickness, the struc- 
ture near interface is more and more affected by the sur- 
face. We conduct simulation on the system with 3 ML 
of Pt and 10 ML of AgCo. Co islands appear only in 
the very early stage, and the final structure is a complete 
layer structure without any Co island. 

Our analysis suggests there exists a critical value for 
the periodic strain of substrate to induce nanostructure 
in thin film. A previous study [p] indicates that the 
imposed strain provides a pre-existing preference wave- 
length for the spinodal decomposition in initial stage, re- 
gardless of how small the strain is. The simulation results 
in this section shows that such preference wavelength 
may be lost in later stage and self-organized nanostruc- 
turcmay fail to realize if the strain is less than a critical 
value. 



E. Film growth: layer by layer 

The above research simulates the annealing process of 
the disordered AgCo phase where atomic diffusion occurs 



in the whole thin film (bulk diffusion). In experiments, 
surface diffusion is much faster than bulk diffusion. In 
the epitaxial growth of thin film, the system is produced 
in layer by layer. If one controls growth conditions care- 
fully to prohibit interlayer diffusion while fully permit- 
ting intralayer diffusion, the resulting nanostructure can 
be markedly different (an interesting example in experi- 
ment is the formation of metastable CuCo alloy in thin 
filmQ). We have simulated such a growth process by 
confining the diffusion to within each layer. The result of 
a system with 10 ML of AgCo and CA g = 0.75 is plotted 
in Fig. 12. No layer structure is observed since interlayer 
diffusion is prohibited. The patterned nanostructure near 
interface penetrates through the whole film to the sur- 
face. This would result in the phenomena of growth on 
seed layer, in which case an as-grown AgCo layer influ- 
enced by the substrate strain and shape in nano-pattern 
acts as seed layer for the forthcoming growth. Since 
AgCo is a system with spinodal decomposition, Ag atom 
is attracted to Ag atom while Co to Co, so that the nano- 
pattern of the as-grown layer is transferred as film grows. 



IV. CONCLUSIONS 

Based on the above simulation results, it is con- 
cluded that patterned nanostructure can be obtained in 
AgCo/Pt/MgO thin film. The main feature of AgCo 
spinodal decomposition in this system is the competi- 
tion between the surface-directed layer structure and 
the strain-induced patterned structure. Strain effect is 
weaker than the surface effect in the current cases. The 
patterned Ag-Co structure occurs only within a short 
range near the interface (no more than 10 monolayers), 
and the thickness of Pt phase should be small enough to 
effectively transfer the strain of dislocation in Pt/MgO 
interface to the AgCo phase to produce nanostructure. 
However, if spinodal decomposition of AgCo phase oc- 
curs within monolayers by careful film growth control- 
ling, surface-directed layer structure will be prohibited 
and the patterned structure will form in a wide range. 
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